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Abstract 

For a particle channeled in the bent crystal planes (axes), the 
phenomenon of "bending dechanneling" is well known, which is 
a particle transition to random state due to centrifugal force. 
We consider an analytical theory of the reverse phenomenon, i.e., 
feeding from a random to channeled state due to centrifugal force 
in a crystal with variable curvature. 

1 Introduction 

For any trajectory of a particle in a crystal a time-reversed trajectory is pos- 
sible. The starting point of a trajectory becomes the final one, and vice versa. 
This leads to the idea of reversibility of the dechanneling processes []J • In the 
depth of a crystal, besides the particles leaving the channeling mode (dechan- 
neling or feeding out) , there may be particles entering the channeling mode 
(respectively, feeding in, also known as volume capture). The mechanisms, 
responsible for these two opposite processes, are essentially the same. Two 
mechanisms changing the particle state in a crystal, namely the scattering 
and the centrifugal effects, are known (see references in Ref. B]). 



"Electronic address: biryukov@mx.ihep.su 



1 



Physics Letters A 205, 343-348 (1995) 



2 



In a bent crystal (or a straight crystal with lattice dislocations), the 
particle state may be changed by a change of a curvature (centrifugal effects). 
The existence of the trajectories starting in the region of channeled states and 
ending outside (or vice versa) can be shown by solving the equation of particle 
motion in a crystal with variable curvature (see, e.g., Ref. |§). Along these 
trajectories the particles leave or enter the channeling mode. The bending 
dechanneling is a well established phenomenon. Its theoretical description 
within the framework of the continuum model 0, is in excellent agreement 
with experiment ||. The same theory inevitably leads to the existence of 
the corresponding feed-in process |7|, for the crystal channels with variable 
curvature (crystal bending or dislocations). 

We shall consider the centrifugal mechanism of feeding-in by the analyt- 
ical methods. The central problem is to find the probability of feeding in a 
general case of the charged particle interaction with the potential of the crys- 
tallographic planes of variable curvature. We shall give below two derivations 
for the probability of feeding: (1) from the consideration of particle dynamics, 
and (2) from the consideration of reversibility of trajectories. 

2 Derivation from dynamics 

Let us consider an efficiency of feeding into a planar channel with increas- 
ing bending radius R; the interplanar potential is U(x). A particle (with 
momentum p and velocity v ) moves as if it were in the effective interplanar 
potential 

U eff (x) = U(x) + ^-x (1) 
with the transverse energy 

E T = ^ + U eff (x) (2) 

where 9 is the particle angle relative to the planes (we denote with 9 C and 
x c the critical angle and transverse position, respectively). A change of cur- 
vature k = pv/R determines, through the relations ([l],|D a corresponding 
change of the transverse energy Et of a particle in an effective potential 
well. For rapid variations of k, when the effective well changes significantly 
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over one period of oscillation A, the term Et looses its meaning, and one 
should solve the dynamical equation, for a given function k(z). 

In practice, A is quite small (~30 fim at 100 GeV), so the crystal curvature 
cannot considerably change along this length (it is not true for dislocations). 
Therefore, the case of a slowly varying curvature (when the "transverse en- 
ergy" concept remains useful) is general enough; it covers all the bent-crystal 
experiments done thus far. In this case, the particle continues to oscillate in 
the bent channel, with all the parameters changing slowly. 

As a particle moves in the crystal, the local curvature k(z) varies with the 
crystal depth z. Instead of the consideration of a variable critical transverse 
energy E c (z), it is more convenient mathematically to keep E c constant and 
to introduce an equivalent variation of Et for any particle. With a small 
instant variation 5k, the transverse energy (TJ) of the particle (in a "new" 
effective well at z+5z) varies, with respect to the critical value E c =U e ff(—x c ), 
by (x + x c )5k. Since both 5k and x are dependent on z, the total increase 
is 5Et = j(x + x c )5k, where x may depend on the fc(2;)-function for rapid 
variation of k(z). For slow variations, one may rewrite it as 

5E T = k' J(x + x c )dz = k'((x) + x c )5z , (3) 

where we have averaged x over oscillation, and kept k' =dk / dz=coiast along 
it. 

Figure |1] helps to understand the bending capture. The effective inter- 
planar potential is shown at some coordinate z and at z + A. A random 
particle may be captured when it is passing over the potential well with a 
transverse energy being in the range from Uq to Uq + AU. The potential well 
gap AU = pvdp/R is due to centrifugal force (d p is the interplanar spacing). 
During oscillation over the well with Ex — Uq, the Et variation relative to 
E c is 

5E T = k\(x) +x c )\ . (4) 

Note that (x)<0 for positive pv/R. If R increases, 5Et < and some 
particles (those with Et + 5Et < E c ) are captured. 

Figure |I| shows three kinds of particles. Particles with Et < Uq + 5Et 
are captured; those with E T > Uq + AU — 5'Et (where 5'Et is variation 
over A/2) are transmitted and interact with the next wall; other particles are 
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reflected by the channel wall. For slow variation the number of transmitted 
particles may be neglected; then we normalize the captured fraction to all 
the particles incident in AU range. For the flat distribution of transverse 
energies in the interval U$ to Uq + AU, the fraction of captured particles is 

5E T R'Xx c f.,{x)\ 

where R' = dR/dz is the curvature radius gradient. For slightly-bent crystals 
the factor in the parentheses is close to unit. Then one can simplify to 

WB ~ lit (6) 
In the harmonic approximation Eq. (|5]) takes the form: 

R'Xx c ( R c \ 

R c being the critical radius. 

In the widely-applied 3-point bending scheme, the curvature is a linear 
function of z: 1 / R(z)=coiast±aoz . Interestingly, this leads to a linear depen- 
dence of wb on R (like in the case of 'scattering feeding' ||, |I0| ) for not too 
strong bending: 

w B = a XR^ (l + &-)*^R. (8) 



x. 



c 



wb is of order of 2X/Lb, where Lb is the length of the bent part of a crys- 
tal. The probability grows with energy, and should be quite sizable, e.g. in 



conditions of the H8 experiment [DJ at CERN, where (§) constitutes several 
per cent. 



3 Derivation from reversibility 

The capture probability can be derived with no use of approximations, like 
the existence of E?, or a slow varying curvature (although these are well 
justified in all practical cases). 

Let us consider a beam with a uniform angular distribution, 1/2$, in- 
cident on a crystal with a curvature gradually increasing along the beam 
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path from to 1/ R c . The fraction of beam channeled will decrease from the 
starting value given by 

/=^t|, 0) 



d p 4 $ 

with the crystal depth, according to: 

f(z) = ^^A B (pv/R(z)) (10) 

where A B (pv / R(z)) describes the reduction of the bent-crystal acceptance 
with depth z. The factor 7r/4 is exact in a harmonic case; for the realistic 
potential one should replace 7r/4 with a slightly different factor ^0.8. 

The number of particles dechanneled over the length 5z, from z to z + dz, 
is equal to 

df(z) _ 2x c n 9 c dA B (pv/R) pv dR(z) 
dz d p 4 $ d(pv/R) R 2 (z) dz 1 ' 

All the dechanneled particles are distributed in the angular range from to 
the bending angle of the crystal. The particles dechanneled over dz are ex- 
iting in the angular range dO = dz/R(z). Therefore the angular distribution 
downstream of the crystal is 

d f - f>( \ d f( z ) _ 2x c ir6 c dA B (pv/R) pv dR{z) 

d6 ~ [Z) dz dp 4 $ d(pv/R) R(z) dz [ } 

where z and 9 are related through the expression d9 = dz/R(z). For an 
abrupt (R'—oo) or very fast variation of curvature, Eq. (|T2| ) overestimates 
the phase density, of course, as we don't take into account a small extra 
spreading, ±9 C , of the dechanneled particles. The microscopic phase density 
of the particles is to be conserved; df/d8 cannot exceed the initial value. 

Let us now consider the same beam incident on the same crystal in the 
reverse direction. That is, the crystal curvature is gradually decreasing along 
the beam path. Now the particles with the upstream parameters (xi,8i) 
equal to the downstream parameters (xf,9f) of the dechanneled particles 
from the preceding case, are captured along the same (reversed) trajectories. 
By consideration, the number of particles which have experienced transitions 
from the channeled states in the former case, is equal to the number of 
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transitions to the channeled states in the latter case (as the trajectories are 
the same). 

Therefore the number of particles captured from the interval d9 is given 
by Eq. (|T2p. Normalizing it to the number of particles incident on the crystal 
in this angular range, d9/2<&, we obtain the capture probability: 

,„ r 9 \ _ Qjfr#(*) _ ™cfi dA B (pv/R) pv dR 

w B {z) - 2$ — - —9 C d(pv/R) - ■ — (13) 

The above quantity is a function of z, through R(z) and R'(z). The derivation 
used just the general arguments, relating the capture probability to the well- 
known function Ab{pv/R) of bending dechanneling. We cannot show the 
equivalence of ( |L3|) and @ for an arbitrary U(x). For a harmonic U(x), 
A B = (1 - R c /R) 2 |]. Then flU) transforms into 

2irx B 6 B A R c \ R C R' x Cx f A R c \ R' 

which exactly equals the result of previous section, Eq. (0); the following 
relations were used here: \=2nx c /9 c , R c = pvx c /2E c = x c /9 2 . The consis- 
tency found between the above two very different derivations, Sects. || and 
H makes us more confident in the theoretical predictions. 

It may be interesting to say that the same reversibility approach has 
been applied to the scattering mechanism, leading to a derivation of the cor- 
responding feed- in probability: ws = tvR9 c /2Ld {Lb is the 1/e dechanneling 
length due to scattering) |T0 . 



4 Integral trapping efficiency 

Crystal may trap particles, with probability (^,|l3]), in the angular range 

A9= dz/R{z) , (15) 

J Z\ 

where Z\ and Z2 are defined by the interval with R' > and R > R c . The 
capture angle A# may be arbitrarily large (with respective loss in efficiency). 
For maximal efficiency, A9 should be within the divergence of the incident 
beam. From the preceding consideration it is already clear that the number 
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of particles fed-in from a divergent beam through the considered mechanism 
equals, in the optimal case, to the number trapped from the same beam at 
the entry face. It is interesting, however, to show explicitly for the harmonic 
potential, from Eq. (|7]) that the integral efficiency for a divergent beam coin- 
cides with the entry-face capture efficiency for the same beam, independent 
of the behaviour of function R(z). 

The total number of particles, captured into the channeling mode from 
the beam with an angular distribution P(9), equals to: 

N B = J w B P(9)d9, (16) 

where one can integrate over z: d9 = dz/R(z). For the uniform distribution 
P{9) = 1/2$, covering the angular region of capture A6 < 2$, we have 

J 1 dzR'Xx c f R c \ Xx c R fdR f R c \ 
B = J^~R~Rd7{ ll) = 2&LJl¥{ ^"J' ( " ,] 

which is independent of the R(z) form. For the potential well deformation 
from Ri = R c to R 2 = R, the number of captured particles equals 

Xx c f R c \ 2 n9 c 2x c f R c \ 2 

which is exactly equal to the efficiency for the entry-face capture. 

This conclusion, derived analytically in the harmonic approximation, is 
true also in the general case of an arbitrary potential, as a consequence of the 
Liouville's theorem. In the absense of dissipative processes the microscopic 
phase density of the captured beam equals the phase density of the incident 
beam. 

The angular acceptance of a crystal deflector could be made much wider 
than the divergence of an incident beam, A# 3> 2$. In this case the intensity 
decreases according to (|P7D, where R\ and R 2 are the bending radii at the 
points, where the beam envelopes are tangent to the crystal planes. 

5 Experimental evidence 



The first experimental indication of the role of the gradient effect in volume 
capture was made at 70 GeV JTJ]. The measurements gave for wb the values 
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factor of 2 higher than those predicted by Eq. (||). The dependence of wb on 
R was measured to be that predicted by Eq. ([5]) (linear in case of Ref. ||12|). 

6 Additivity of the two mechanisms of feed- 
in 

In bent crystals both mechanisms, bending and scattering, may contribute 
to the feeding-in. In the following we present a qualitative analysis, which 
indicates that the two mechanisms contribute in an additive way. 

The 'scattering feeding' is efficient only at the tangential points. The 
particles pass over the potential well with the initial transverse energies being 
in the range from Uq to Uq + AU. Having scattered during this oscillation, 
the particles get some distribution }\Et) with a (interesting to us) tail to 
smaller < Uq. Those particles which have E? < E c , are trapped in the 
stable channeled states. Their number is 

N s = f(E T )dEr, (19) 
Jo 

The curvature variability modifies E c \ then the trapped fraction is 

N = / f(E T )dE T = N s + / f{E T )dE T (20) 

Jo JE C 

Since normally 5E C <C E c over one oscillation A, we have 

N = N S + f(E c )5E c = N S + N B . (21) 

The term Nb = f(E c )5E c ps 5E c /AU is the fraction of particles trapped 
through the bending mechanism, as it is defined in Section 

For positive 5E C (the well grows), the two contributions are summed: 
iV = Ns + Nb- For negative 5E C (the well shrinks), the contributions are 
subtracted: N = Ns — Nb- 

7 Conclusion 

Although the processes considered cannot increase the beam deflection ef- 
ficiency, they contribute to a wide angular acceptance 9 C ) crystal de- 
flector. The deflection efficiency (always smaller than in the regular case of 
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the entry-face capture) of such deflector can be designed, and be varied in a 
broad range of values. Moreover, this efficiency is independent of the incident 
beam divergence (in contrast to the regular case). This may be valuable, for 
instance, for the proposed [JT3] measurements of CP-violation with use of Kg 
and Kl meson beams at CERN, where beam attenuation, stability, and low 
background are important issues. 

Obviously, the secondary particles produced in collisions with crystal nu- 
clei can be captured into the channeling mode, from the crystal bulk, through 
the feed- in processes only. This may be applicable, e.g., in some ideas of crys- 



tal usage for the experiments in particle physics [14 



As the energy increases, the rate of transitions from the random to the 
channeled states changes. For a given crystal, the rate ws of the scattering- 
induced transitions decreases as l/(pv) 3 ^ 2 since the scattering vanishes [[K] . 
In contrast, the rate % of bending capture, Eq. (||), grows as (pv) 1 ^ 2 since 
the centrifugal effects strengthen. The ratio of the two feeding mechanisms 
changes in favour of bending with the energy increase. 

At sufficiently high energies (above ~1 TeV), the processes of feeding-out 
and feeding-in in bent crystals are predominantly defined by the centrifugal 
effects. These processes are described by the continuum model, which is a 
well established analytical theory. Furthermore, these centrifugal processes 
are easily controlled, through the proper design of the i?(z)-function. 
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Figure 1: Schematic picture of bending capture. The effective interplanar 
potential U e ff(x) is shown at some coordinate z (o), and at z + A (•). There 
are three kinds of particles: (1) captured, (2) reflected, (3) transmitted and 
interacted with the next wall. 



